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We present the results of measurements of the thermal conductivity of the quasi one-dimensional 
spin S—1/2 chain compound SrCu02 in the temperature range between 0.4 and 300 K along the 
directions parallel and perpendicular to the chains. An anomalously enhanced thermal conductivity 
is observed along the chains. The analysis of the present data and a comparison with analogous 
recent results for Sr2Cu03 and other similar materials demonstrates that this behavior is generic 
for cuprates with copper-oxygen chains and strong intrachain interactions. The observed anomalies 
are attributed to the one-dimensional energy transport by spin excitations (spinous) , limited by the 
interaction between spin and lattice excitations. The energy transport along the spin chains has a 
non-diffusive character, in agreement with theoretical predictions for integrable models. 
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I. INTRODUCTION 



The physics of one-dimcnsional (ID) magnetic systems 
has been of interest for some time but particularly signif- 
icant progress, both in experiment and theory, has been 
made during the last several years. The related research 
activities have been growing because the quantum na- 
ture of low-dimensional low-spin systems promises a rich 
variety of phenomena to be explored. Among other fea- 
tures, the transport of energy in ID magnetic systems is 
expected to be highly unusual. A number of models de- 
scribing one-dimensional systems are integrable, imply- 
ing, for instance, conservation of energy current and, as a 
consequence, ballistic energy propagation and divergent 
thermal conductivity. The question whether the energy 
transport is diffusive jOiJjallistic is currently under active 
discussion for atomim a as well as spincTO ID systems. 

The ID Heisenberg 5=1/2 systems with antiferromag- 
netic coupling between adjacent spins are of particular 
interest. It is well established that elementary excita- 
tions in such systems are not S=l magnons but S—1/2 
quantum solitonsJ3 commonly called spinous. Spinous 
can only be created or destroyed in pairs but they may 
be treatedr-igs free particles, at least at large spatial 
separation.Ercl The interaction of spinous, as quantum 
solitons, with structural defects and other quasiparticles, 
such as phonons, is poorly investigated both theoreti- 
cally and experimentally, although in classical Heisenberg 
and Ising ID magnetic systems, the interaction of soli- 
tons with defects, phonons, andjapaanons demonstrates 
a number of interesting featurestjllil (for a review, see 
Ref. |l|). 

In order to probe the features of energy transport 
in 8=1/2 Heisenberg chains, only a limited number of 
experiiaents, such as inelastic neutron scatteringcB or 
NMR,tZl is available. In case of a sufficiently strong cou- 



pling between the phonon and the spin subsystems, the 
energy transport mediated by magnetic excitations can 
also be monitored via experiments probing the thermal 
conductivity. Although the spin-lattice interaction as 
well as the influence of defects and interchain interac- 
tions make the spin system non-integrable, an anoma- 
lous behavior of the spin-mediated thermal conductivity 
may, nevertheless, be expected. Establishing the heat 
transport carried by spin excitations may, in addition, 
reveal information concerning magnetic defects and the 
spin-lattice interaction. The spin-lattice interaction is 
very important in ID magnetic system, because it leads, 
to modifications of the spectrum of spin excitationsEs 
and, under special circumstances, to the formation of new 
phases, such as the spin-Peierls dimerised state. Besides 
defects and phonons, interchain interactions are also an- 
ticipated to influence the energy transport in the spin sys- 
tem and hence may be probed by thermal-conductivity 
experiments. 

In previous work concerning the thermal conductiv- 
ity k(T) in quasi-lD S=l/2 Heisenberg systems, obser- 
vations of some spin- mediated heat transport, in addi- 
tion to the dominant-phonon contribution, have.. been 
reported for KCuFa,^ CuGe03,E] and Yb4As3.El The 
most pronounced effects, however, have been observed 
in (Sr,Ca)i4Cu2404i, containing doubl e . kg . Cu-0 lad- 
ders with 180° Cu-O-Cu bond anglcs.BEl The ther- 
mal conductivity along the ladder direction exhibits an 
anomalous double-peak temperature dependence.|Jt has 
fairly well been established by a detailed analysiscj that 
this two-peak feature in n(T) of (La,Sr,Ca)i4Cu2404i is 
caused by two main contributions, one due to phonons, 
responsible for the low-temperature peak, and the other 
due to itinerant spin excitations in the ladders, causing 
the additional high-temperature maximum. 

In order to test the above mentioned conjectures con- 
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cerning energy transport by spin excitations in 5'= 1/2 
ID Heisenberg systems, we have investigated the low- 
temperature thermal conductivities of two materials con- 
taining 180° Cu-0 chains, nav aehf SrCu02 and Sr2Cu03. 
The crystallographic unit cellsEjE^ of the two compounds 
are schematically drawn in Fig. N_ Mic intrachain ex- 
change coupling J of 2100-3000 KEHea in these two ma- 
terials is of the same order of magnitude as the interac. 
tions along the ladder legs in Sri4Cu2404i ( J=1510 K).E2l 
The compounds considered here exhibit different ar- 
rangements of the same type of chains. The structure 
of Sr2Cu03 contains Cu-0 chains with a verii-.small in- 
terchain interaction J' {a = J' / J ~ 10~^).l3 SrCu02 
is built by double Cu-0 chains forming Cu-0 ribbons 
containing Cu-Cu zig-zag chains (see Fig. 0). The inter- 
action between the two chains occurs via 90° Cu-O-Cu 
bonds providing a weak (|a| = 0.1-0.2) ferromagnetic 
interaction.cil For comparison, the chains forming the 
double-leg ladders in (Sr,Ca)i4Cu2404i are connected via 
180° Cu-O-Cu bonds ("mngs"), therefore for the ladders 
a = 0.55 is rather largeEj The spin excitation spectrum 
for the single-chain compound Sr2Cu03 is gapless. Be- 
cause the interaction between the chains forming the Cu- 
O ribbons in the SrCu02 is-fcrromagnetic, the excitation 
spectrum remains gapless,EHl in contrast to the gapped 
spectrum of spin ladders where the interaction between 
the chains is antiferromagnetic. 

In this paper, we present new experimental results on 
the thermal conductivity of SrCu02 along the main crys- 
tallographic directions in the temperature region between 
0.4 and 300 K. These data are analyzed together with, 
those for Sr2Cu03 measured and presented previously.c^ 
We offer a detailed analysis of phonon and spinon con- 
tributions to the thermal conductivity along different 
crystallographic directions. We show that the thermal 
transport perpendicular to the chain direction is pre- 
dominantly phononic at all temperatures, whereas along 
the chain direction the spinon contribution is dominant 
at elevated temperatures. The scattering of phonons is 
mainly via other phonons, at dislocations, and at sam- 
ple boundaries. The spinon thermal conductivity may 
reliably be evaluated only in the high temperature re- 
gion above 60 K where it is predominantly limited by 
the spin-lattice interaction. At very low temperatures, 
the spin-lattice interaction appears to be rather weak 
and therefore the intrinsic spinon thermal conductivity 
cannot be observed. 

For clarity we introduce an auxiliary notation of crys- 
tallographic directions. For both compounds we denote 
the direction along the chains as c' (c for SrCu02 and 
h for Sr2Cu03), the direction perpendicular to the CuO 
squares as h' (a for SrCu02 and c for Sr2Cu03), and the 
direction in the CuO square plane but perpendicular to 
the chains as a' (see Fig. p. This notation is compatible 
with the crystallographic axes of the spin-ladder system 
(Sr,Ca)i4Cu2404i, for which a', 6', and c' correspond to 
the a, b, and c directions, respectively. 



SrgCuOg SrCuOs 




SrgCuOg SrCuOg Sri4Cu2404i 




FIG. 1. The schematic unit cells of the crystal structures 
of SrCu02 and Sr2Cu03. For each material, the layout of the 
Cu-O chains is emphasized. For comparison, the structure of 
the ladder in Sri4Cu2404i is also shown. 



II. EXPERIMENTAL DETAILS 

The single crystals of SrCu02 and Sr2CuOs— .were 
grown by a traveling-solvent floating-zone method.E3 The 
high structural quality of the crystals was confirmed by 
neutron diffraction experimentsBJ 

From each crystal, three specimens in the form of rect- 
angular bars with typical dimensions 2.5x1x1 mm'^ and 
with the longest dimension parallel to each principal axis 
were cut with a thin-blade diamond saw. The thermal 
conductivity was measured by a standard steady-state 
heat-flow technique, fixing the sample at one end to a 
heat sink. The temperature gradient along the sample 
was produced by a heater (1 kil chip resistor) glued to 
the opposite end of the sample, and monitored in over- 
lapping temperature ranges by a matched pair of Ru02 
thermometers at T < 5 K and by differential Chromcl- 
Au + 0.07% Fe thermocouples at higher temperatures. 
The temperature difference between the thermometers 
was typically 1.5% of the absolute mean temperature of 
the sample. The estimated accuracy of ± 10% of the ab- 
solute value of k{T) is mainly caused by the uncertainty 
of the sample geometry, the relative error being of the 
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order of 0.5% of the measured values of k. 



III. EXPERIMENTAL RESULTS 

The resuhs of the thermal conductivity measurements 
for SrCu02 3lc£; presented in Fig. ^, together with our pre- 
vious resulta23 on the thermal conductivity of Sr2Cu03. 
For both materials k{T) reveals similar general features. 
In particular, the thermal conductivities along the direc- 
tions perpendicular to the chains (^a', i^b') are similar not 
only in their features of the temperature dependence, but 
even in absolute values. At the same time, the thermal 
conductivity along the chains is distinctly higher than 
k(T) perpendicular to the chains in both cases. This 
is particularly evident at temperatures above the k(T) 
maxima. At room temperature, this excess conductivity 
for the double-chain material SrCu02 exceeds that for 
the single-chain compound Sr2Cu03 by almost a factor 
of two. In the temperature region between 1 and 6 K, 
where magneticj)lia^e transitions have been reported for 
both materials J^3~E3 no anomalous features in k,{T) have 
been observed. Since these transitions |arp-jndicated by 
anomalies of the magnetic specific heat,E30 the interac- 
tion between the spin and the lattice systems as well as 
the spin contribution to the measured thermal conduc- 
tivity must be negligible in this temperature region for 
both materials. 

Because both SrCu02 and Sr2Cu03 are electrical insu- 
lators, electrons do not participate in the heat transport. 
Below, we thus analyze the experimental data assuming 
that the main heat carriers are phonons and spin excita- 
tions. 



IV. ANALYSIS OF THE PHONON 
CONTRIBUTION 

A. The model and fitting procedure 

In general, the behavior of the phonon thermal con- 
ductivity Kph(T) is determined by scattering of phonons 
by structural defects, phonons and other quasiparticles. 
At temperatures T > Qjj, where Qd is the Debye 
temperature, the phonon-phonon Umklapp-processes (U- 
processes) usually dominate and it often turns out that 
Kph oc T~", where a ~ 1. The mean free path of the 
phonons is very small and comparable to interatomic dis- 
tances. At temperatures below the probability of U- 
processes decreases exponentially and the phonon mean 
free path increases accordingly. The maximum mean free 
path of phonons, normally reached only at very low tem- 
peratures, is limited by the size of a sample. At this, 
so called Casimir limit, the thermal conductivity may be 
calculated by Kph = l/SCphvLh, where Cph is the lat- 
tice specific heat, v is the average sound velocity along 
the chosen direction, and Lf, is a constant related to the 
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FIG. 2. Temperature dependences of the thermal conduc- 
tivities of SrCu02 and Sr2Cu03 along the a', b', and c' axes. 
The dotted lines represent estimated limits of the thermal 
conductivity due to the finite size of the samples. The solid 
and dashed lines represent different evaluations of the phonon 
contribution to described in the text. 



sample size. Hence, for isotropic materials, Kph oc 
at low temperatures. These two sources for scattering 
of phonons are always present in any material, generat- 
ing some typical features of K;ph(T) including a distinct 
maximum at Tmax < 6_d/10. Various defects, such as 
impurities, dislocations, stacking faults etc., also partici- 
pate in phonon scattering, thereby reducing the thermal 
conductivity, especially in the region of Tmax- Various 
quasiparticle excitations such as magnons, electrons, ex- 
citons etc. also influence the thermal conductivity by 
scattering phonons and by providing additional channels 
of heat transport. 

In order to identify the dominant mechanisms of 
phonon scattering and to estimate their relative impor- 
tance, we have fitted the experimental data by using the 
Debye approximation of the phonon spectrum combined 
with the relaxation-time approximation for calculating 
the thermal conductivity. The model assumes the same 
group velocities and relaxation rates for phonons of dif- 
ferent polarizations. A more accurate and more compli- 
cated analysis, taking into account possible differences in 
the scattering rates for phonons of different polarisations, 
would require a detailed knowledge of elastic constants, 
which is still missing. The phonon thermal conductivity 
is thus calculated as 



where lo is the frequency of a phonon, t{uj,T) is the corre- 
sponding relaxation time, is the Debye temperature, 
and X — huj/k^T. As a further simplification we assume 
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that all mechanisms of phonon scattering act indepen- 
dently. In this case, 



(2) 



where each term corresponds to an individual inde- 
pendent scattering mechanism. For instance, 



v/Lb 



(3) 



for phonon scattering by sample boundaries, 

= (4) 



for phonon scattering by point defects (Rayleigh scatter- 
ing), and 
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for phonon-phonon U-processes. As indicated in Eq. (g), 
there may be other contributions to r~^, such as phonon- 
phonon N-processes, phonon-spin scattering, etc. In the 
first step only the three main mechanisms (boundaries, 
point defects and U-processes) were considered for all 
samples. In those cases where the fitting of the data 
with a reasonable accuracy failed, additional terms rep- 
resenting relevant scattering processes were introduced 
in Eq. (|). 

For the fitting procedure we employed the Levenberg- 
Marquardt algorithm. The quantities L;,, A, B, and 6, 
defined in Eqs. (^) to (||) were treated as freely adjustable 
parameters. The criterion of selecting additional terms 
in Eq. (0) was the minimization of the mean-square devi- 
ation We used the @d values of 441 K for Sr2Cu03 
(Ref. g) and 357 K for SrCu02 (Ref. |3|). The mean 
sound velocity v was calculated from the values of 0d, 
using the equation v = Qjj (ks/fi) {6n^n)~^^^ , where n 
is the number density of atoms. 



B. Phonon thermal conductivity perpendicular to 
the chain direction 

As may be seen in Fig. |2[ the thermal conductivities 
of SrCu02 and Sr2Cu03 in the directions perpendicular 
to the chains exhibit the temperature dependence that is 
typical for phonon heat transport. This includes k oc 
at high temperatures and approaching k oc at low 
temperatures, the latter being expected if the phonon 
mean free path Zph is constant. Since Iph turns out to 
be shorter than the smallest sample dimensions even at 
0.4K (for SrCu02), the influence of defects cannot be 
neglected. 

It turns out that the three mechanisms initially in- 
cluded in Eq. (||) were not enough to fit the data in the 
entire covered temperature range. However, after con- 
sidering various other possible scattering processes, very 
good agreement between the experimental and the calcu- 
lated Ka' (T) and Hb' (T) curves in the temperature region 
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T (K) 

FIG. 3. A comparison of the experimental (open circles 
and diamonds) and calculated (solid lines) Mi{T) along the di- 
rections perpendicular to the chains. For clarity, the curves 
for different samples are scaled in order to achieve a reason- 
able separation between them. 



between 0.4 and 100 K was achieved by including two 
additional terms in Eq. (||). The first one. 
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where C is a fit parameter, is usually attributed to 
phonon scattering by strain fields of dislocations .e3 This 
term proved to be essential for fitting the temperature 
dependence of the thermal conductivity at T < Tmax- 
The second term is of resonant typeC^ 



= D- 



, .2^2 



{l-cg^{LJo,T)) 



(7) 



where uq is the resonance frequency, c is the fractional 
concentration of scatterers and 13 is a factor characteriz- 
ing the strength of the resonant scattering. The function 
g{ujo,T) = (p_ — p+), describing the difference in popu- 
lation of the upper (p_|_) and lower (p_) states was taken 
as g{LjjQ,T) = tanh(?i(jJo/2fcBr), corresponding to a two- 
level system with singlet upper and lower levels. The 
parameters D, c, and loq are again fit parameters. A res- 
onant term of this type was initialiy. proposed for phonon 
scattering by magnetic impurities ,cd but was also success- 
fully used to explain the scattering of acoustic phonons 
by flat optical phonon modes.eS 

The fit values of the adjustable parameters Lj,, A, B, 
C, D, loq, c, and b are listed in Table |[ The k(T) curves 
calculated using these values are shown in Fig. ^, together 
with the experimental data. 

For a sample with a rectangular cross-section 5'cs, a 
fully diffusive scattering of phonons at the sample bound- 
aries leads to Lb = The expected values of Lb 
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compatible with the dimensions of our samples are of the 
order of 1 mm, i.e., close to the fit values given in Table |. 
The small values of Lb may possibly be caused by large 
defects, such as microcracks present in the samples. 

The fitting procedure gave values of the parameter A 
in Eq. (|) equal to or less than 10"''^ s^. With A that 
small, the influence of point defects on the phonon ther- 
mal conductivity is negligible in comparison with other 
processes of phonon scattering. 

For materials with v atoms in the unit cell, the pa- 
rameter b in Eq. (j^) is expected to be given by & ~ 21^^/^ 
(Ref. ^3|). The fit values of b presented in Table | are 
in reasonable agreement with this expectation. The pa- 
rameter B, characterizing the strength of U-processes, 
depends in a complex manner on elastic and anharmonic 
constants of a material. 

The parametefi C, associated with dislocation scatter- 
ing, is given byE^ 

C = fcnd6V, (8) 

where rid is the number of dislocations per unit area, b 
is the Burgers vector, 7 is the Griineisen constant, and k 
is a constant of the order of 10~^. Assuming that j ^ 2 
and b ~ 4A, we get ~ 10^'' m~^. 

The origin of a resonant scattering of phonons in 
SrCu02 and Sr2Cu03 is not a priori clear. However, 
the fact that c defined in Eq. (Q) turns out to be close 
to 1 means that the number of scattering centers is of 
the order of the number of unit cells in the sample, 
which severely limits the choice of scattering processes 
and excludes, for example, dilute magnetic impurities as 
scatterers. One possible origin of such scattering could 
be the interaction of acoustic phonons with a flat op- 
tical phonon branch, the situation that was promoted 
in Ref. ^ Another possible source of resonant scat- 
tering are excitations of segments of spin chains of fi- 
nite lengths, created by the inclusion of nonmagnetic 
defects. In this latter case, the expression for g{u!o,T) 
in Eq. (|^) should be changed to describe the degener- 
acy of the ground and excited states of a chain seg- 
ment. For chain segments with an even number of spins, 
which are exoected to prevail in S = 1/2 AFM Heisen- 
berg chains,C3 the low-energy resonance corresponds to 
a singlet-triplet transition and, therefore, glujQjT) = 
{1 — exp{huJo/kBT))/{l + 3exp{hujo/kBT)). It appears, 
however, that choosing between g{ujQ,T) valid for either 
a singlet-singlet or a singlet-triplet excitation has only a 
minor effect on the resulting quality of the fit and the 
values of the parameters in Eq. (|^). Therefore, at this 
point, we cannot distinguish these two possibilities. 

As one can see in Fig. |3[ the experimental data slightly 
deviate from the approximation above 100 K. The magni- 
tude of the deviations increases with increasing temper- 
ature and, at room temperatures, it is as large as 1.5-2 
W (m K)~^. Experimentally, heat losses by radiation 
result in measured values of k(T) higher than the intrin- 
sic thermal conductivity. This effect is often observed at 



high temperatures for samples with poor thermal conduc- 
tivity, if no special precautions are taken. Estimates have 
shown, however, that in our experimental setup the heat 
losses via radiation are at least an order of magnitude 
too small to explain the excess conductivity and hence 
the observed deviation should be attributed to an intrin- 
sic mechanism. Since with increasing temperature the 
heat is mainly transported by high-frequency phonons 
for which the Debye model is not a good approxima- 
tion, the deviations may in principle be attributed to 
shortcomings of the employed simple model. Another 
source of high-temperature deviations may be an ad- 
ditional heat transport by optical phonons. We recall 
that high-temperature deviations of the same magnitude 
can also be discerned in our KadT) data for the spin- 
ladder system (Sr,Ca)i4Cu2404i.c3 The similarity of the 
high-temperature behavior of K.a'.b'{T) of spin-chain and 
spin-ladder materials seems intriguing and may hint to a 
common origin of the feature such discussed later, 

diffusive energy transport via the spin system. 

C. Phonon thermal conductivity parallel to the 
chain direction 

In contrast to the successful interpretation of the ther- 
mal conductivity perpendicular to the chains by consid- 
ering phonon heat transport alone, we failed completely 
in obtaining a reasonably good fit of the anomalous T- 
dependence of Hq' at temperatures above T,„ax with a 
similar approach. In order to emphasize the anoma- 
lous feature in Kc' (T) , we show the anisotropy ratio 
(3 — He' / Hb' in Fig. 0. For comparison, (3 = I i^a' ra- 
tios for the spin-ladder system Sri4_j:Caa;Cu2404i (a;=0, 
2) are included in Fig. ^. At low temperatures the 
anisotropy is relatively weak and varies little with in- 
creasing temperature. A distinct rise of /3 is observed 
above 30-40 K. 

The obvious similarity of the three systems contain- 
ing similar Cu-0 chains with different interchain inter- 
actions suggests a common cause for this anomaly. The 
anomaly might be due to drastic changes in the phonon 
heat transport above Tmax, but we argue below that 
this is not the case. Indeed, in the vast literature on 
properties of Sr2Cu03 and SrCu02 no observation of 
any phase transition above 6 K has been reported and, 
therefore, no drastic changes in elastic interatomic in- 
teractions at T > Tmax is expected. Without such a 
change affecting the lattice, the chain-like structure of 
the lattice itself cannot explain the strong temperature- 
dependent anisotropy of the thermal conductivity above 
Tmax- To illustrate this point, we show in the inset of 
Fig. 1^ the anisotropy of the phonon thermal conductivi- 
ties for IIg2Cl2 and Te (data from Refs. ^ and ^ respec- 
tively), two chain- type materials with pronounced elas- 
tic anisotropics. IIg2Cl2 and Te were chosen as examples 
because in both cases K(r) varies as at high tem- 
peratures, i.e., the phonon-phonon scattering dominates 



5 





6p 




^ 4 - 


10 










^ 2- 








1 . 

1 


00. 





SrigCa2CU2404T" 


, ,,| 


■ Sri4CU2404i— ' 


1 ' 
I > 
1 1 


SrCuOg /"^ 


/ 

/ Sr^CuOa, 


/ 1 
~ / Lj 


^-^Te _ 




HgzCij \^ 
■r:■■:■■:^ - 



Sri4CU2404i 
Sr-| 2C32CU2404-| 




10 

T (K) 

FIG. 4. Temperature dependences of the anisotropy 
ratios k^i / for SrCu02 and Sr2Cu03, and n^i / hi^i for 
(Sr,Ca)i4Cu2404i (data of Ref. [j^). The inset illustrates the 
relative change of the anisotropies with temperature above 
Tmax for the same compounds and also for Hg2Cl2 (data of 
Ref. ^) and Te (data of Ref. 



above Tmax- This demonstrates that even for chain-type 
materials with strong elastic anisotropies, the anisotropy 
of the phonon thermal conductivity above T^ax is only 
weakly temperature-dependent. 

From the point of view of elastic properties, the spin- 
chain cuprates are rather ordinary, weakly anisotropic 
materials. For example, the anisotropy of the sound ve- 
locities of Sri4Cu2404i is found to be very weak.cZl There- 
fore, as we have already argued in Ref. ^ well above 
Tmax one may expect not only a weak temperature de- 
pendence of /3, but also a rather small anisotropy of the 
phonon-phonon scattering, and therefore /3 ~ 1. 

On the other hand, the high-temperature features of 
Kc' {T) do correlate with features in magnetic ppiperties. 
In Fig. H we reproduce the results of Zhai et aLcS for the 
dynamic magnetic susceptibilities Xci'^iT) of SrCu02, 
Sr2Cu03, and Sri4_a;Caa;Cu2404i (x — 0, 2) in the c'- 
direction, measured at w = 10 GHz. By comparison of 
Figs. ^ and ||, the correlation between the thermal con- 
ductivity and the rria*;netic susceptibility along the chain 
direction is obviouala and indicates that the anomalous 
Kc' (T) is of magnetic origin. 

Based on the arguments presented above, the increase 
of P{T) for SrCu02 and Sr2Cu03 with increasing temper- 
ature is attributed to an additional, non-phononic chan- 
nel of heat transport along the chain direction. Most 
naturally, this excess thermal conductivity is associated 
with an energy transport via spin excitations of the spin 
chains. This additional heat transport is large enough to 
be separated reliably from the phonon contribution. 

In order to achieve this separation of the spin- mediated 
contribution to the thermal conductivity from the ex- 
perimental data on k^', some assumptions about the 
phononJpackground have to be made. As we have argued 
before]^ the phonon-phonon scattering is most likely al- 
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FIG. 5. Imaginary part of the dynamic magnetic sus- 
ceptibility x"(10GHz,T) (reproduced from Ref. ksl), for the 
compounds that are discussed here. 



most isotropic. For the evaluation of the relaxation rate 
via phonon-phonon scattering along the chain direction, 
one can therefore use the average values of the param- 
eters B and 6, obtained by fitting k{T) along the other 
directions (see Table |). With this assumption, we have 
fitted the data for k.^' (T) at low temperatures (T < 10 K) 
where the influence of spin-mediated heat transport, as 
argued above, is negligible. The fitting procedure was 
the same as for Ka' and k;,/ but now, only the parameters 
L, A, C, and D were kept free. For the parameters c 
(concentration of resonant scatterers) and ujq (resonance 
frequency), again the average values obtained by fitting 
Ka' (T) and Kb' (T) were used. The best fits are shown in 
Fig. y as solid lines and the corresponding fit parameters 
are listed in Tabic |. 



D. Evaluation of the heat transport in the spin 
system 

The spin contribution to the heat transport along 
the chain direction was extracted by subtracting the cal- 
culated phonon contribution from the experimental data 
of Kc'(T). The resulting Ks{T) dependences are shown 
in Fig. H as open symbols. It is these data that are 
used for the subsequent analysis of spinon heat transport. 
However, the way of evaluating Ks, described above, 
strongly relies on the assumption of isotropic phonon- 
phonon scattering. Alternatively, one may consider that 
the anisotropy of the phonon thermal conductivity does 
not completely vanish for T > Tmax, but is similar in 
magnitude as for the case of chain materials Hg2l2 and 
Te discussed above (see the insert of Fig. The limiting 
case is then that the anisotropy of the phonon thermal 
conductivity is constant and equal to that in the region 
of Tmax- This situation is illustrated in Fig. |^ where the 
dashed line represents the phonon thermal conductivity 
along the c'-axis, averaged along the a' and b' -axes and 
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FIG. 6. Spinon thermal conductivity for SrCu02 and 
Sr2Cu03. The solid lines are calculated using Eqs. ( |l^ ) and 
(^|). The shaded areas demonstrate possible errors caused by 
the uncertainty of the estimate of the phonon thermal conduc- 
tivity. The inset shows the ratios between the effective and 
the intrinsic spinon thermal conductivities, calculated using 
Eq. M. 



scaled to the value of Kc'(71nax)- The result of subtract- 
ing this calculated phonon contribution from the experi- 
mental Kc' (T) values fixes the lower boundary of possible 
values of Kg (T) and coincides with the lower edge of the 
shaded areas in Fig. I These areas indicate the possi- 
ble uncertainty of Ks{T), caused by the ambiguity of our 
evaluation of the phonon parts of the thermal conductiv- 
ities. The uncertainty is large at low temperatures, but 
at higher temperatures it is small enough to allow for a 
quantitative analysis of the thermal conductivity via spin 
excitations. 



to be nearly independent of temperature. Using, for the 
specific heat, the low-tejaperature (T ^ J/ks) result of 
a Bethe ansatz solutiontJ for a uniform chain of N spins 
S = 1/2 



Cs 



^Nk% 
3 J ' 



the thermal conductivity in the diffusive regime is 



21,2 



a k 



(9) 



(10) 



Here, Ug is the number of spins per unit volume. Eq. ( p^ ) 
predicts rather small values for (of the order of 
1 W (m K)"-'^ at 300 K) decreasing linearly with decreas- 
ing temperatures. The observed Ks{T) along the c'-axis 
is much higher than predicted by Eq. ( p^ and, at least 
at T > 80 K, it increases with decreasing temperature. 
This is a clear sign of non-diffusive energy transport via 
spin excitations along the chains, in agreement with the- 
oretical predictions. It is worth noting that also the spin 
diffusion constant Ds has experimentally been found to 
be strongly enhanced in Sr2Cu03.E3 This observation was 
considered as indicating the ballistic nature of spin trans- 
port. 

In contrast to the c'-direction, the high-temperature 
excess thermal conductivity along the a' and b' direc- 
tions, discussed in section IV B, is in fair agreement with 
Eq. (|l^ both in absolute values and with respect to 
its temperature dependence. We therefore argue that 
the possible energy transport perpendicular to the Cu- 
O chains via the spin system is small and diffusive. On 
the contrary, the energy transport via spins along the 
chains is substantial and relies on the ballistic propaga- 
tion of spin excitations, interacting with lattice imperfec- 
tions and other quasiparticles. 



V. ANALYSIS OF THE SPIN CONTRIBUTION 

A. Heat transport via diffusion 

An important question which can be answered imme- 
diately is whether the observed energy transport via spin 
excitations is diffusive or ballistic. For uniform ID spin 
systems the energy transport is expectei-to be diffusive 
for spin S >l/2 and balhstic for 5=1/2.E£| However, the 
diffusive energy transport in 3—1/2 chains may be re- 
stored by, for example, interchain interactions or intra- 
chain next-nearest-neighbor interactions^ In case of dif- 
fusive energy transport, the thermal conductivity is given 
by Ks = DeCsi where Cs is the spin specific heat per unit 
volume and D_e is the energy diffusion constant. In the 
high-temperature limit, De for a quantum spin system 
is given by De = KJ[S{S + l)]^/^a^/?i, where a is the 
distance between neighboring spins and if is a constant 
which|4n different theoretical models, varies between 0.25 
and 3B In the low-temperature region, De is expected 



B. Spinon heat transport along the chain direction 

As mentioned in the introduction, the elementary exci- 
tations of a uniform Heisenberg ID AFM S" = 1/2 system 
carry a spin S — l/2n Therefore the low-temperature 
thermodynamic behavior of such an ensemble can be ac- 
counted for by using the usual Fermi-Dirac distributiop. 
/ — {exp{e/kBT) + 1)^^ with zero chemical potential,El 
where e is the spinon energy. The dispersion relation for 
spinous is given byi3 

Jtt 

e(fc) = — |sinfca| . (11) 

For spinous, the possible values of the wave vector k are 
restricted to one half of the Brillouin zone. The low- 
temperature specific heat, calculated using the general 
equation 
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FIG. 7. The specific heat of the S=l/2 AFM Heisenberg 
chain, calculated using Eq. ( |l^ (solid line) and the Bethe 
ansatz solution (dashed line, data from Ref. pT]) . 



agrees with the low-temperature result of Eq. (||) . As il- 
lustrated in Fig. 0, the high-temperature specific heat cal- 
culated with Eq. (|l^) qualitatively reproduces the main 
high-temperatupe,Jeatures of the results of the Bethe 
ansatz solution.E3E3 The disagreement between the two 
models is not important for our analysis since our re- 
sults have been obtained in the low temperature region 
fceT/J < 0.15. 

Applying the Boltzmann-type approximation, the 
thermal conductivity of a spinon system is 

■K/2a 

2nsa f df 
K-s = j -Q^evslsdk, (13) 



where Vs = h~^de/dk is the group velocity and Ig is 
the mean free path of a spinon. Cast in another but 
equivalent form, 

jTv/2kBT 

^s^^T J ^^Ue,T)d., (14) 



where x — e/fceT. If several independent scattering 
mechanisms have to be considered, the spinon mean 
free path can be represented as l~^{e,T) = ^ l~l{e,T), 
where each ls^i(e,T) term corresponds to an indepen- 
dent scattering channel. This is analogous to the case 
of phonons, considered in Eq. (^. 

The most straightforward way to analyze our Ks{T) 
data would be analogous to that employed for the phonon 
conductivity, i.e., by choosing a set of ls,i{s,T) terms in 
the equation for the mean free path of spinous and fitting 
the data of Ks{T) to Eq. (p^. However, the energy and 
temperature dependences for most of the plausible scat- 
tering mechanisms have not yet been worked out. That 
is why we treated the Kg (T) data by inserting the spinon 
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FIG. 8. The temperature dependences of the spinon 
mean free paths for Sr2Cu03 and SrCu02, calculated using 
Eq. (0. The solid lines are fits using Eq. (|l|). The shaded 
areas have the same meaning as in Fig. ^. 



mean free path IsiT), averaged over all e, at a given tem- 
perature. In that case, ls{T) may be taken out of the 
integral in Eqs. (|l^) and (|lj) and thus may be calcu- 
lated from the Ks{T) data, using the known values of J 
for the materials investigated here. For these calcula- 
tions, J/ks = 2100 K (Ref. ||) and 2620 K (Ref. |3|) 
were used for SrCu02 and Sr2Cu03, respectively. The 
resulting values of ls{T) for SrCu02 and Sr2Cu03 are 
shown in Fig. ||. 

Here again, the shaded areas reveal the possible range 
of Is values caused by the uncertainties in evaluating the 
phonon thermal conductivities. In Ref. it was noted 
that IsiT) for Sr2Cu03 may be approximated by 

Is' = Isp + l;J = ^,Texp(-T7T) + Lj, (15) 

where the two terms were attributed to spinon scatter- 
ing by phonons and defects, respectively. Here, the pa- 
rameter Asp characterizes the strength of spin-lattice in- 
teraction, T* is related to the minimum energy needed 
to produce a single spinon-phonon Umklapp-process, 
and Lsd is the mean distance between magnetic defects 
effectively scattering spin excitations in a spin chain. 
It may be seen that with Eq. (|l5|) and the appropri- 
ate parameters, ls(T) may be well approximated not 
only for Sr2Cu03 but also for SrCu02. The values of 
Asp=7.1xl0^ m-iK-i, r*=177 K, L^d=8.1xl0-* m 
for Sr2Cu03, and ^,^=6.7x10^ m^^R-i, r*=204 K, 
is(i=2.42x 10^^ m for SrCu02, were obtained by fitting 
the data of IsiT) with Eq. (|^. The fits are shown as 
solid lines in Fig. ^. The parameter values for Sr2Cu03 
are slightly different from those presented in Ref. |3^, 
because of a somewhat different way of evaluating the 
phonon thermal conductivity. The agreement between 
the experimental data and the fit is very good above ap- 
proximately 50 K, but at lower temperatures the uncer- 
tainty in extracting IsiT) becomes far too large to draw 



any reliable conclusion concerning its temperature de- 
pendence and the involved scattering mechanisms. 

The fit values of the parameters A^p and T* which 
are related to the spin-lattice interaction arc similar in 
magnitude for the two compounds, in spite of the differ- 
ent chain arrangements. The values of T* scale with the 
Debye temperatures, supporting the idea that the term 



Ijp in Eq. ( |l5[) is indeed related to the spin-lattice in- 
teraction. It appears that the main difference between 
the two materials is in the defect scattering, which is 
three times more efficient for Sr2Cu03 than for SrCu02. 
The low-temperature values of Ig for Sr2Cu03 have been 
founded to be in fair agreement with the mean distance 
between neighboring bond-defects causing local changes 
in the magnetic cmipling between Cu spins, as extracted 
from NMR data.Ea Unfortunately, similar NMR results 
revealing the distance between bond-defects are absent 
for SrCu02. Our data suggest an about three times 
smaller concentration of these defects in SrCu02 than 
in Sr2Cu03. 



C. Influence of spinon-phonon scattering 

For the two investigated materials, the evaluation of 
the spinon part of the thermal conductivity has been 
achieved only in a limited temperature region, charac- 
terized by 0.015 < T/J < 0.15. Two mechanisms 
of spinon scattering seem to be effective, namely the 
temperature-dependent spinon-phonon interaction and 
the temperature-independent scattering of spinous by de- 
fects. The significance of other scattering mechanisms, 
if any, is negligible. For example, normal processes of 
spinon-spinon scattering cannot change the total mo- 
mentum of the spinon system and, therefore, should not 
contribute to the thermal resistance. The spinon-spinon 
Umklapp-processes are expected to play a serious role 
only at temperatures T ~ J/fce . 

Supposing that the validity of Eq. extends beyond 
the region 0.015 < T/J < 0.15, we may estimate the ex- 
pected spinon thermal conductivity at lower and higher 
temperatures using Eq. (O) . The result is shown as solid 
lines in Fig. ^ Since I s{T) is constant and equal to Lgd 
at very low temperatures, Ks{T) is expected to vary lin- 
early with T as Ks = nsak^TrLsdT /2ih. The phonon ther- 
mal conductivity varies as at low temperatures and 
hence the spinon contribution should progressively dom- 
inate with decreasing temperature. This expectation is 
obviously not met by our low-temperature results. First 
of all, no linear contribution to the total measured ther- 
mal conductivity has been identified and, second, no dis- 
tinct feature in k{T) is observed at the Neel temperature 
Tjv, indicating a negligible spinon contribution k^. This 
suggests that at some temperature below 30 K the spinon 
contribution deviates from the predicted behavior shown 
in Fig. H and by 6 K, spinous are excluded from heat 
transport. At this point, a discussion of the spin-lattice 



interaction is in order. 

Although the spin-lattice interaction reduces both the 
phonon- and the spin-related thermal conductivities be- 
cause of scattering processes involving both types of 
quasiparticles, some degree of it is needed for the spin- 
related heat conduction to. be observable in a thermal 
conductivity experiment The energy provided by a 
heater generates only phonons, and the spin-phonon in- 
teraction is needed for an energy transfer from the lattice 
to the spin system. The effective spinon thermal conduc- 
tivity Ks.qff] which is accessible in a coupled spin-lattice 
system, isE3 



Kg tanh(ylLsample /2) 



where Lsampic is the sample length. 



1/2 



A = 



(16) 



(17) 



and Tsp is the spin-lattice relaxation time. For very short 
spin- lattice relaxation times (j'^p'^ — ^ ^50)5 '^s,eff — f^s- 

In 

the opposite case of a very weak spin-lattice interaction 
{t^p^ — !■ 0), however, off = and the thermal trans- 
port by spinous cannot be observed in an experiment, 
regardless of how large the intrinsic Ks of the sample is. 

If the first term in Eq. ( p^ is indeed dictated by the 
spin-lattice interaction, it is possible to establish the tem- 
perature region in which the observed Ks,cff equals the 
intrinsic k^. For this purpose, the ratio Ks^cbJj^s was 
calculated by assuming Tap = hp/vs (see Eq. ( p^ ) ) and 
Vs — Jair/ih, valid for low-energy spinous which domi- 
nate at T ^ J/kB- The result is shown in the inset of 
Fig. ^. Although the extrapolation of Eq. (p^ ) to lower 
temperatures may not exactly be correct, it is neverthe- 
less clear (see the inset of Fig. |^) that at very low tem- 
peratures the spin-lattice interaction turns out to be too 
weak for the spinon thermal conductivity to be observed. 

From our discussion it is obvious that the interaction 
between the spin and the phonon system should show 
up both in spinon and phonon heat transport. How- 
ever, other mechanisms of phonon scattering, mainly the 
phonon-phonon interaction, may be stronger and mask 
the scattering of phonons by spin excitations. This may 
be the reason why a term corresponding to the phonon- 
spinon scattering does not have to be included in Eq. (|^) 
for a reasonable description of Kp\i{T). On the other 
hand, a phonon-spinon relaxation rate fps(T), averaged 
over the entire phonon spectrum at a given temperature, 
may roughly be approximated by 



' ps 



ph 



(18) 



where at a fixed temperature T, Nph{T) and Ns{T) are 
the average numbers of phonons and spinous per unit 
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volume, respectively. At high temperatures, where spin- 
phonon interaction is important, Eq. ( |l8|) for the phonon- 
spinon scattering rate gives oc T exp{~T* /T), i.e., 
the same temperature dependence as is characteristic for 
the phonon-phonon interaction. That is why, even if the 
influence of phonon scattering by spin excitations is not 
negligible, it is difficult to separate it from the phonon- 
phonon interaction. 

In the present discussion, we do not consider the pos- 
sibility of a mutual drag between spin and lattice exci- 
tations, which may be considerable under certain condi- 
tions, as described in Ref. 



VI. SUMMARY 

In this paper the thermal conductivities of the spin- 
chain compounds SrCu02 and Sr2Cu03 were studied. 
Although the crystal structures of the two compounds 
are different in the sense that the former contains linear 
and the latter zig-zag Cu-0 chains, the thermal conduc- 
tivity of both materials is remarkably similar. In particu- 
lar, the heat transport in the directions perpendicular to 
the chains is dominated by phonons, but along the chain 
direction and at high temperatures there is a substantial 
excess contribution related to the transport of energy by 
spinous. 

The phonon thermal conductivity is analyzed employ- 
ing a Debye-type approximation. The main sources of 
phonon scattering are phonons at high temperatures and 
lattice defects, presumably dislocations, at low temper- 
atures. The spin-phonon interaction is not seen in the 
phonon heat transport, most likely because it is masked 
by other scattering processes. 

An eventual energy transport via the spin system in the 
directions perpendicular to the spin chains is found to be 
small and, if any at all, of diffusive character. On the 
contrary, the spin-mediated energy transport along the 
chain direction is large and propagating ballistically. The 
contribution of spin excitations (spinous) is reliably sin- 
gled out along the chain direction at high temperatures 
where it exceeds the phonon contribution quite substan- 
tially. Our analysis based on a fcrmion model suggests 
that the expected infinite spinon thermal conductivity is 
limited by the influence of defects and phonons. 
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TABLE I. Parameters of the fitting of k{T) data to Eq. (0). 



Parameter 


SrCu02 


SrCu02 


SraCuOs 


SraCuOa 


SrCu02 


SraCuOa 




a'-axis 


6' -axis 


a'-axis 


&'-axis 


c'-axis 


c'-axis 


L (10-^ m) 


0.32±0.1 


0.25±0.02 


0.97±0.5 


8±6 


0.63±0.04 


2.5±2 


A (10-*2 s3) 


< io-*5 


< io-« 


< 10-^5 


< 10-^5 


< io-« 


< 10-45 


B (10-^** s K-^) 


11.1±0.5 


13.5±0.5 


7.6±0.5 


9.0±0.5 


12.3 


8.30 


b 


2.7±0.2 


2.8±0.1 


3.5±0.3 


3.3±0.2 


2.75 


3.4 


c iw~^) 


7.2±1 


5.8±0.3 


9.4±1 


10.6±0.7 


3.9±0.2 


5.5±0.2 


D (10'' s-i) 


4.6±2 


2.7±1 


5.6±3 


6.5±2.5 


3.5±1 


7.7±0.2 


^0 (10^' s-i) 


lO.Oil 


8.3±1.5 


10.2±1 


9.6±1 


9.2 


9.9 


c 


0.86±0.03 


0.90±0.01 


0.76±0.03 


0.87±0.03 


0.88 


0.815 


Temperature region 


T <100 K 


T <100 K 


T <100 K 


T <100 K 


T <12 K 


T <13 K 
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